Superradiant terahertz emission by dipolaritons 
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Dipolaritons are mixed light-matter quasiparticles formed in double quantum wells embedded 
in microcavities. Due to resonant coupling between direct and indirect excitons via electronic 
tunnelling, dipolaritons possess large dipole moments. Resonant excitation of the cavity mode 
by a short pulse of light induces oscillations of the indirect exciton density with a characteristic 
frequency of Rabi flopping. This results in oscillations of classical Hertz dipoles array which generate 
supperradiant emission on a terahertz (THz) frequency. Resulting THz signal may be enhanced using 
the supplementary THz cavity in the weak coupling regime. 



Introduction. Generation and frequency modulation 
of terahertz radiation is one of major technological chal- 
lenges nowadays Applications of terahertz sources 
span from communication technologies to medicine and 
security. The existing THz emitters are based on a large 
variety of physical principles Q. One possible solution 
is a conventional solid state oscillator based on a high 
frequency Gunn or tunnelling diode The operating 
frequency of such oscillator is limited to the microwave 
region and the lower boundary of terahertz range. Sec- 
ond, quantum cascade lasers (QCL) allow for coherent 
emission of radiation in a terahertz range 0, i]. Ex- 
ploiting the multiple photon emission from intersubband 
transition of wide quantum wells Q, QCL cover the up- 
per boundary of THz range with relatively high power of 
emission (up to 50 mW). QCL operate at cryogenic tem- 
peratures which strongly limits their application area. 
The laser driven terahertz emitters form a third group 
of THz sources, where femtosecond optical pulse illumi- 
nating the semiconductor structure leads to oscillation 
of carrier density which generates terahertz radiation 0- 
Q. Similar effects can be observed in crystals with a 
large second-order susceptibility where time- varying po- 
larization induces THz emission (lo| . Finally, the wide 
group of THz emitters are free electron based sources like 
klystron or free electron laser. These are powerful but 
quite bulky sources. Optimization of spectral character- 
istics, size, efficiency, cost and operation temperature of 
terahertz sources is one of the priority objectives of mod- 
ern optoelectronics. Here we propose a compact optical 
to terahertz radiation convertor allowing for efficient fre- 
quency modulation and operating at high temperatures. 

Theoretical proposals for THz sources in solid state 
and semiconductor physics present a rich diversity. One 
example is the use of carbon nanostructures, in partic- 
ular, carbon nanotubes (CNT's) [HI, UH, or Aharonov- 



Bohm quantum rings j 1 31 ] - Another area where a possi- 



bility of THz generation has been widely studied theo- 
retically is polaritonics — the interdisciplinary research 



area at the boundary of solid state physics and quan- 
tum optics 0, [HJ]. Polaritonic devices are based on 
the strong coupling between excitons in semiconductor 
quantum wells and confined photons. Proposals include 
the polariton based THz emitters where signal is gener- 
ated from transitions between upper and lower polariton 
branches [H, HH, EH and transition between 2p and Is ex- 
citon states, the latter one strongly coupled to the cavity 
mode Q . Recent theoretical studies also include the pro- 
posal of bosonic cascade lasers realized due to the multi- 
ple THz photon emission by excitons (exciton-polaritons) 
confined in a parabolic potential trap (20l |. Our present 
proposal is based on the recent experimental realization 
of dipolaritons — cavity exciton-polaritons characterized 
by large dipole moments. 

The mechanism for terahertz signal generation de- 
scribed in this Letter relies on the beats between spa- 
tially direct and indirect excitons. Spatially indirect ex- 
citons are composed of electrons and holes situated in 
separate quantum wells (QWs) (2l| . Coupled QW sys- 
tems where separation of electrons and holes and energy 
splitting between direct and indirect excitons can be con- 
trolled by applied electric fields have been studied in the 
recent decades fiiWzl ] . The radiative lifetimes of indirect 
excitons are usually much longer than those of direct ex- 
citons due to the lower electron-hole overlap. Another 
important feature of indirect excitons is their large dipole 
moments in the normal to QW plane direction resulting 
in strong exciton-exciton interactions |25| . 

Recently, it has been shown that exciton polaritons 
and spatially indirect excitons can be intermixed in the 
biased semiconductor microcavities with embedded cou- 
pled QWs @, S, Hfl • In these structures new quasiparti- 
cles being linear superposition of cavity photon (C), di- 
rect exciton (DX) and indirect exciton (IX) appear. They 
form three exciton-polariton modes, namely, the upper 
polariton (UP), middle polariton (MP) and lower polari- 
ton (LP) modes. These modes may be characterized by 
large dipole moments which is why they are referred to 
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as dipolaritons |8(. 

Here we study theoretically the effect of beats between 
dipolariton modes due to the tunnel coupling between 
direct and indirect excitons. The combined tunnelling 
and cavity-exciton coupling effects result in Rabi flopping 
between two dipolariton modes and lead to harmonic os- 
cillations of IX density. Since indirect excitons repre- 
sent elementary dipoles, these density oscillations result 
in emission of THz frequency radiation whose charac- 
teristics can be tuned by the applied bias and pumping 
conditions. This mechanism of THz radiation is similar 
to downshift optical-to-THz frequency conversion used in 
laser driven terahertz emitters and allows for high power 
generation. The essential difference of the generation 
scheme we propose from the previously studied schemes 
is in the use of the superradiance effect appearing due 
to the fact that Rabi oscillations in the dipolariton sys- 
tem are coherent. Rabi oscillations can be sustained by 
the cavity for several tens of picoseconds. This strongly 
improves the quantum efficiency of the emitter. 

The model. The structure we consider represents two 
QWs separated by a barrier which is sufficiently thin to 
allow for resonant electron tunnelling (Fig. [T]). The elec- 
tron wave function is shared between two QWs in this 
case. The optical microcavity is tuned to the wider QW 
exciton resonance (the left QW in Fig. [T]), while the right 
QW remains decoupled from the optical pump. Pulsed 
pumping creates electron-hole pairs which form direct 
excitons. The bias applied in growth direction induces 
mixing of direct and indirect exciton states. Coupling 
of the cavity photons to direct excitons resulting in ap- 
pearance of two polariton modes has been extensively 
studied @, [l5| . The coupling and anticrossing of DX-IX 
resonances tuned by electric field is also a well-known 
effect in GaAs/AlGaAs QW structures [H,^. Both ef- 
fects combined result in the appearance of dipolaritons 
studied recently by Cristofolini et al. Q . 

The correct treatment of a real dipolariton system 
should in principle involve both coherent and decoher- 
ent parts for the Hamiltonian, H = H C oh + Pi dec- The 
Hamiltonian corresponding to the coherent part reads 

iicoh =fiu)ca^a + hwoxb^h + hjj IX 6^6 + — (a)b + Pa) 

(1) 

- ^-(b ] c + $b) + Ptf + P*a, 

where a' , a, b\ b and c\ 6 are creation and annihilation 
operators for cavity photons, direct excitons and indirect 
excitons, respectively. Here Swc, fu^DX and hwix denote 
cavity mode, direct exciton and indirect exciton energies, 
and first three terms of Hamiltonian describe energy of 
bare modes. The fourth and fifth terms correspond to 
coupling between modes, where HQ (Rabi splitting) [32| 
denotes the coupling constant between photon and direct 
exciton, and the tunnelling rate corresponding to DX-IX 
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FIG. 1: (Color online) Sketch of the system, (a) Double 
quantum well (DQW) heterostructure with resonantly cou- 
pled electron levels and a hole in the left quantum well (QW). 
The right QW has larger bandgap than left QW tuned by the 
presence of indium alloy or width of the well, (b) Electron- 
hole bilayer with schematic picture of coupled spatially indi- 
rect and direct exciton which form dipolariton. 



coupling is hJ. The last two terms in Hamiltonian cor- 
respond to the optical pumping of cavity mode with rate 
P(t) = Po(t)e tu}pt , where top is a pumping frequency. 
In the following considerations we omit summation over 
space indices since we are mostly interested in temporal 
dynamics of the system. The decoherent part of Hamilto- 
nian Wdec in general can contain the terms corresponding 
to different dissipative processes, including finite lifetimes 
of the modes and pure dccohercncc. The dynamics of the 
system can be calculated from Master equation for the 
density matrix, 



= [Hcoh,p] + tp, 



(2) 



where C denotes a Linblad superoperator responsible for 
accounting of incoherent processes 0, 0. Its detailed 
form is given in Supplementary material, section A [3lj . 

In the following Letter we arc interested in the large oc- 
cupation numbers of the modes, and mean-field approxi- 
mation is eligible. The dynamics of occupation numbers 
of the modes can be found from the Master equation ^ , 
where operators are treated as classical fields [the details 
of this derivation are given the Supplementary Material, 
section A] 
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where (...) = Tr{...p} denotes averaging using the den- 
sity matrix p, and occupation numbers for modes can be 
found as nc = (a)(a)* , n DX = (b)(b)* &ndn IX = (c)(6)*. 
Here we applied the rotating wave approximation (a,-) — > 
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FIG. 2: (Color online) Dynamics of the dipolariton system 
subjected to pulsed optical pumping which shows oscillations 
of indirect exciton (red line, nix) and direct exciton (blue 
line, riDx) density. The coupling constants are equal to hJ = 
hQ, = 6 meV. (a) Oscillations of density in regime I, where 
cavity mode is detuned from direct exciton resonance by Sn = 
— 10 meV. Pumping frequency is chosen as hA = 11.5 meV, 
and we plot oscillations for electric field F = 0.95-Fb, where 
Fq is field corresponding to IX-DX resonance (e.g. Fq — 12.5 
kV/cm in Ref. [g). The green area schematically represents 
optical pulse with duration At = 1 ps and pulse edge defined 
by 0.1 ps transient time (pump, scaled intensity). The inset 
shows the long-term term oscillations of nix(t) and nox(t) 
density which are in antiphase. The period of oscillations is 
equal to T = 0.65 ps. (b) Oscillations corresponding to regime 
II, where detuning is equal to 8n — 3 meV. For the same 
electrical field and pump intensity (HA = 0) the magnitude 
of signal is higher comparing to regime I. Inset shows highly 
damped long-term oscillations in the second regime. 
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and introduced the photonic 5c 



ujc ~ udx 



and excitonic 5j = u>ix — Wnx dctunings. The pumping 
rate reads as P(t) = Po(t)e~' lAt , where A = uip — ujc 
denotes detuning of the coherent optical pumping. The 
damping rates of the modes are defined by parameters 
7, = 2ir/ri, i = C, DX, IX. Typical lifetimes of the 
modes are tq ~ 3 ps, tdx ~ 1 n s and tjx ~ 100 ns. 

Additionally, we derived dynamic equations for occu- 
pation numbers for the modes going to higher order of 
mean- field theory and verified validity of the system (|34l) 



(|36p for chosen pumping conditions. The former model 
represents the general description of dipolariton system 
and allows for accounting of various decoherent processes 
including pure dephasing [see Supplementary informa- 
tion, section B]. 

Results and discussion. We calculate numerically the 
dynamics of the system with coupled quasiparticles sub- 
jected to a picosecond pulsed optical pumping. The pres- 
ence of mixing terms between different modes implies the 
oscillating behavior similar to Rabi flopping in the clas- 
sical model of a two-level system subjected to a time- 
varying field. It is important to note that our system has 
several characteristic frequencies. They are governed by 
the exciton-photon coupling strength f2 and detuning Sq, 
the IX-DX coupling strengths J and detuning Sj. Addi- 
tionally, the pumping frequency cup governs the efficiency 
of the pump. Varying these characteristic frequencies one 
can control the frequency, magnitude and damping rate 
of indirect exciton density oscillations. 

While the coupling constants and J are dependent 
on the geometry of the structure and can hardly be tuned 
for a given sample, the detunings between modes Sj and 
5q are strongly sensitive to the applied electric field and 
the incidence angle of the cavity pump 0] . Tuning these 
parameters one can bring the oscillating dipole system to 
different regimes. If the cavity mode is far-detuned from 
IX-DX anti-crossing the light-exciton coupling is weak, 
which will be referred to as the regime /. If the detun- 
ing Sq is small the strong intermixing of IX, DX, and C 
modes takes place, which corresponds to the regime II. 

The behavior of the system in the regime / is shown 
in Fig. Ufa) for the detuning 5^ = — 10 meV. We observe 
antiphase oscillations of IX and DX densities with de- 
caying amplitudes. The inset in Fig. [2ja) shows a zoom 
of long-standing phase-locked oscillations which last for 
several tens of picoseconds. For the electric field F ap- 
proximately equal to the resonant field of IX-DX modes 
anti-crossing Fq the frequency of oscillations is given by 
v rj J/2tt = 1.45 THz. The time dependence of indirect 
exciton density oscillation nix (t) , calculated numerically 
from system (f3"4")) - (|3"fj|) , can approximately be fitted with 
analytical function as 



nix(t) = n? x cos 2 (wi/2)e" t/T 



(6) 



where ui = 2ttv is the frequency of oscillations, n° IX de- 
notes the magnitude of oscillations which decreases in 
time with the damping rate r _1 . Tuning the electric 
field, the frequency of generated oscillations changes in 
the range of several THz due to its dependence on IX- 
DX detuning (Fig. [5]). Other important parameters of 
the system such as the amplitude of oscillations and di- 
mensionless oscillation quality factor £ defined as ratio of 
magnitude to the decrement of oscillations can be tuned 
by variation of pumping conditions, as well as by the ap- 
plied electric field [see Supplemental Material, section D, 
for details]. 
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FIG. 3: (Color online) Frequency of indirect exciton density 
oscillations as a function of dimensionless electric field cal- 
culated for oscillations in regime I (solid line) and regime II 
(dashed line). 



Decreasing the detuning between exciton and photon 
modes Sq one can bring the system into the regime II. 
In this regime the IX density oscillations are observed 
as well, while their quality factor is different. Strong in- 
teractions of the cavity mode with the IX-DX resonance 
result in a higher magnitude of nix oscillations than in 
the regime / (Fig. EJb)). The damping rate of these os- 
cillations is also larger. Clearly, the regime of strong cou- 
pling between all modes is advantageous for the pulsed 
pumping regime and it allows for high power generation. 
On the other hand, the regime / is preferential for the 
long-standing signal generation providing higher quality 
factor £, while lower amplitude of oscillations. 

We have shown that due to the coupling between 
modes IX and DX densities oscillate with a THz fre- 
quency. This infers that a dipolariton is an oscillating 
dipole, with a dipole moment in z direction d z changing 
periodically from value d z — do corresponding to an in- 
direct exciton (IX) to d z = corresponding to a direct 
exciton (DX). One can define the total dipole moment 
of the system as D z = Njxd z , where Nix = nix A is 
the number of indirect excitons within the area A illumi- 
nated by the pumping light. Since we have shown before 
that after initial transient regime the density of indirect 
excitons nix(t) is decaying harmonic function of time, 
the total dipole moment of the system can be found as 



D z (t) = d Q ni X Acos 2 (ujt/2)e- t/r . 



(7) 



Here nix is the maximum density of indirect excitons. 

The mechanism of THz radiation by an oscillating den- 
sity of indirect excitons is similar to the photo-Dcmbcr 
effect where oscillations of electron-hole plasma on a sur- 
face of semiconductor generate THz signal [tl, and 
other laser driven terahertz emitters. However, the sys- 
tem of dipolaritons has some important advantages over 
previously studied systems, namely: (1) the better tun- 
ability of the system allowing for fast modulation of the 
THz emission frequency; (2) improved spectral charac- 
teristics of THz signal which can be controlled using ap- 
plied voltage and pumping conditions; (3) possibility to 



achieve a high output power as it is described below. 

The total intensity of the far field radiation emitted by 
a classical Hertz dipole can be found from [34| 



D 2 



67rerjc J 



(8) 



where D z is the total dipole moment of the dipolariton 
array, erj is vacuum permittivity and c is the speed of 
light. For the particular case of an array of harmonic 
dipole oscillations the intensity is 
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(9) 



where do = eL is a dipole moment of indirect exciton 
with L being separation between centers of QWs. Here 
for simplicity we did not consider the damping part e~'/ r 
of oscillations, which can be added straightforwardly. It 
will result into damping of total power of emission, which 
is why to achieve stable cw radiation one therefore needs 
to use sequence of pump pulses. Similarly to the conven- 
tional case of elementary dipole emitter, the polar pat- 
tern is given by Ig ~ sin 2 8, where 9 is an angle between 
direction of radiation and the axis of the structure (see 
Supplementary information, section E, Fig. 3). 

It is important to note that total emitted intensity is 
proportional to the square of the density of indirect ex- 
citons and is dependent on the pump intensity in non- 
linear way. This is a manifestation of the superradiance 
effect: due to the interference of coherent in-phase 
oscillations of elementary dipoles the output power is en- 
hanced. This effect is sensitive the quality factor of the 
cavity: the longer Rabi oscillations persist the stronger 
the amplification effect is. Superradiance is a specific fea- 
ture of the dipolariton THz emitter, which makes it more 
efficient than any existing laser-to-THz convertor. 

For hJ — 6 meV and the typical distance between 
QWs of L = 12 nm, the power of the THz emission of 
one elementary dipole formed by a dipolariton is Iq = 
1.8 x 1CP 17 W = 18 aW. The total power emitted is given 
by this quantity multiplied by the number of elementary 
dipoles Nix- For the typical concentration of indirect 
excitons achievable experimentally of nix = 10 10 cm~ 2 
and a 60 nm. diameter of the pumping spot one obtains 
I to t ~ 1.4 /i W. By embedding a stack of double quantum 
wells nuQw — 4 in each microcavity and using several 
sets of cavities on one chip one can obtain the output 
power similar to one of quantum cascade lasers. 

Furthermore, to improve the efficiency of THz radia- 
tion, the system can be placed in an external cavity tuned 
to the THz frequency. It can be created using metallic 
mirrors or using an inductor capacitor circuit [fl [38| ■ The 
sketch of the system is shown in Supplementary informa- 
tion [section E, Fig. 3]. The efficiency of the emission 
would be further increased by the Purcell factor of the 
external cavity ((J. For the chosen frequency 1.45 THz 
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the experimentally achieved Purcell factor of an inductor- 
capacitor cavity is about Fp = 17 [38 1. 

Conclusions. We have shown that the system of 
dipolaritons can provide an efficient tunable source of 
THz radiation. In this system, due to the direct-indirect 
exciton coupling, the optically excited dipolariton system 
exhibits density oscillations with a subpicosccond period. 
The dipole oscillations lead to the superradiant emission 
of radiation in the THz frequency range. The spectral 
properties as well as the power output of this system are 
expected to be strongly improved with respect to the 
existing laser induced THz emitters. Additionally, we 
propose a way to enhance the radiation efficiency using 
a supplementary THz cavity in the weak coupling regime. 
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Excellence in Polaritonics" and FP7 IRSES project "PO- 
LATER". O. Kyriienko acknowledges the support from 
Eimskip Fund. We thank Sven Hofling and Jeremy 
Baumberg for useful discussions on the subject. 
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Supplemental Material: Superradiant terahertz emission by dipolaritons 



A: Derivation of dynamic equations (lowest order of mean-field theory) 

In the following section we present a derivation of Eqns. (3)-(5) which describe the dynamics of dipolariton system. 
For the macroscopic occupation of the modes one can treat fields as classical and calculate the dynamics of the system 
using mean value of their creation (or annihilation) operators. Using the Hamiltonian (1) and Master equation (2), 
one can write the equation of motion for mean value of operator d as 

hd t (a)(t) = iTr{a[p,U coh ]} = iTr{p[H coh , a}}, (10) 

where we used the cyclic properties of the trace. The commutator in a RHS using explicit form of the Hamiltonian 
reads 



loco) a + uj DX b r b + ljixc* c + ^(d f & + ^d) - — (ftc + tfb) + Pa) + P*d,d 



(11) 



and can be calculated part by part. Naturally, only first, fourth and penultimate terms contribute to the equation. 
Using standard commutation relations for bosons, one gets ujc[a)a, a] = -wpa, §[d^6, a] = — ^b, and P[dt,d] = —P. 
The coherent part of the equations for (b) (t) and (c) (t) can be written in the same way. 

The incoherent processes can be treated using Lindblad superoperator In particular, the accounting of the life 
time of cavity mode can be included to the Eq. (JTUJ) as 

d t (a)(t) = Tr{d^^} = Tr{a^(2apa) - a) ap - pa) a)}. (12) 

The first term gives 7cTr{pd^dd} = 7cTr{pdd^d} — 7cTr{pd}, where we used the relation (d^d)d = d(d^d — 1). 
The second term reads as — ^Trjpdd'l'd}, and third term yields — ^Tr{pd^aa} = — ^Trjpdd^d} + ^Tr{pd}. The 
summation of all terms results into incoherent part of dynamics in the form 

dt(a){t) = —(&)• (13) 

The decoherent part of (b) (t) and (c) (t) can be derived in the same fashion, leading to the system of three coupled 
equations (3)-(5) in the main text which describes dipolariton system. An analogous system of equations can be 
derived from the three coupled oscillators Hamiltonian, which is shown in Section C. 



B: Derivation of dynamic equations (higher order of mean- field theory) 



In the previous section we derived equations which describe dipolariton system in the lowest order of mean-field 
theory. Here we generalize it to the next order and derive equation for the occupation numbers of the modes correlators. 
In the following derivation we will treat coherent and incoherent processes separately. 

Coherent part. We start from the Hamiltonian given by Eq. (1). The time dependence of the occupation numbers 
and correlators can be found using relation 

Ti{a}aid t p{t)}, 



d t rii(t) 

where d^ represents operator of different type. The equation can be rewritten in the form 

d t rii(t) = ^{alailpiHcoh}}, 



(14) 



(15) 



and simplified using the explicit expression for coherent part of the Hamiltonian. To describe the dynamics we need 
to write system of six coupled equations for occupation numbers of the modes nc = (d^d), rinx = (b b), njx = (e'e) 
and correlators a a b = (a^b), on, c = (b^c) and a ac = (a^c). Using mean- field approximation, we derive equation each 
type of occupation number and correlator. For instance, the occupation number of indirect excitons can be written 
as 



d t n IX {t) = iTi{p 



ojca'a + ujuxb^b + ujix^c + — (a)b + b^a) — ~^(b^c + c^b), 



}■ 



(16) 
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The first four terms of commutator are zero, while last term gives 



b^c, ffc 
with final equation 



c^b, cJc 



J - ([S+a, fi f ]c + c f pc,c] + [c t 6,c t ]c + &$b,$[) = -~ (S^-c^) , (17) 



d t nix(t) = JIma 6c , (18) 



where we used relation ImA = (A — A*)/2i. One can write equations for nc and nnx in the same fashion. 

Derivation of dynamic equations for correlators require more efforts, but can be done straightforwardly. For example, 
the equation for correlator a a b = (b'6) reads as 



d t a ab (t) = iTv{p 



LOca^a + uijjxb^b + uiix^c + -^o^b + — + ^^b, a)b 



}• (19) 



Considering terms one by one, the first gives wc[a'«, a^b] = weba^ and second term yields uiDx[b b\ aJb] = —uiDX^b- 
Naturally, third and forth terms give zero, while the fifth term is equal to Q,/2\ab\ a^b] = —£l/2(a^a — b^b). The sixth 
term reads — J/2[wc, a^b] = J/2 a^c, and last term from correlator in Eq. (|19j) is zero. It is convenient to perform 
the same analysis for correlator (5.6^) = a* ab , and combine both expressions. This leads to dynamic equation for the 
real part of the correlator 

d t Rca ab (t) = (u>nx - Lu c )lTs\a a b - ^Ima„ c , (20) 

where relation Re^4 = (A + A*)/2 was used. Equations for dtctbc and dta ac can be derived analogously. 

Finally, assuming the presence of coherent pump, we need to add three equation for (a), (6) and (c) given by Eqns. 
(3)-(5) in the main text of the article. This leads to complete system of Eqns. ([28)) - (j36|) for dynamics of dipolariton 
system written in the end of the section. 

Decoherent part. The Master equation for density matrix is a powerful tool which allows one to account for various 
incoherent processes. It can be introduced in the model using Lindblad superopcrator. In particular, the effect of 
temperature can be modelled using phonon bath with chosen Bose distribution function, as was done in Refs. [l[ 
and Q • In the following considerations we do not consider fc-space dynamics of the system and energy relaxation of 
dipolariton modes. However, we do need to account for dissipation processes coming from leakage of cavity mode out 
of the mirrors, as well as radiative and non-radiative recombination of excitons. 

Dissipation processes can be modelled with Lindblad superopcrator written for any operator of type i in the 
form 

p. 

C-disP = «y (2apa t - a)ap - pa) a) , (21) 

where Ti = tvji represent the decay rate of the mode, i = C, DX, IX. Reminding the Master equation, = 
\H, p] + Cp, one can write the incoherent part of dynamic equation for occupation numbers and correlators. For 
instance, in the case of cavity photon occupancy it reads 

dtneit) = Trja^a = Tr{a^a-^-(2apa^ — a)ap — pa) a)} = 7cTr{p(a^d^aa — a^aa^a)}, (22) 

where to write the latter expression we used the cyclic properties of the trace. Applying the mean-field approximation, 
(nc) = (o)a) = Tr{ficp} and (a)a)aa) = Tr{nc(nc — l)p}> it can be rewritten as 

d t n c (t) = 7 C Tr{p(n c (n c - 1) - h 2 c )} = -^ c nc- (23) 

The same considerations are valid for dtnnxit), dtnixif) and correlators. 

The lifetime of the mode is not the only cause which affects dynamics of the system. Indeed, the coherence in the 
system can be spoiled due to dephasing even for infinitely long-lived atomic systems. There the main effect of loosing 
coherence can be introduced by pure decoherence term, which acts solely on off-diagonal terms or "coherences" , while 
live diagonal terms are unaffected. Usually written for two- level system (atom, qubit etc), the Lindblad operator for 
pure coherence in general form reads Q 

t decP = -iT^ ec )[a z ,[a z ,p]}, (24) 
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where a z = [<j\cf]/2, is a fermionic creation operator, and r( dec ) denotes rate of pure dephasing. In particular, 

for two-level system it thus can be rewritten as CdecP = i r< 2 {pzpGz ~ P) [H- I n our system we want to account for 
decay of coherence of the interaction between different modes which have bosonic statistics. Therefore, we can map 
angular momentum operators into bosonic operators using Holstein-Primakoff transformations a z = (s — d^d), 
where s = 1 is spin of the particle. Then Lindblad operator (|24|1 can be rewritten opening the commutator and reads 

£decP = iT^ dec ' (2a [ apa) a — d^dd^dp — pd^dd^d) . (25) 

Plugging (|25[) into equation for occupation number 

d t n c (t) = 7 (dec) Tr{d f d (2d t dpd t d - dWdp - pdWd)} = (26) 

one can check that it gives zero contribution to the lifetime of the mode. Performing the same analysis for correlator 
(e. g. a ab ) we straightforwardly derive 

dta a b(t) = 7^ ec ' l Tr{d^6 (2d^dpd^d — aJaa^ap — pa)aa)a)} = 7^ ec - ) Tr{pd^S(2dd^d^d— d^dd^d— dd^ddt)} = — 7^ ec ' ) (d^S) . 

(27) 

Thus, coherence between mode is destroyed for large damping rates 7^°, an d index corresponds to the type of 
correlator we are interested in. Similarly, one introduce pure decoherence for a bc and a ac correlators. 
Finally, merging coherent and decoherent parts, we come to the system of equations 



dnc/dt = mma ab - lc n c - 2Im{P* (d) } , (28) 

dn DX /dt = -ttlma ab - Jlma bc - "fDxn D x, (29) 

dnix/dt = JIma bc - -fn IX , (30) 

d{Rca ab ) J 7 ah (dec)^ j f ™m-, , ,x 
-Q t = -o n Ima ab - — Ima ac — Rca a fc - 7^ ; Rea afc - Im{P {b)\, (31) 

d(Rea bc ) il j bc (dec) 

— = djlma bc - — lma ac — Rea bc - 7^ 'Rca hc , (32) 

d{Rca ac ) = _ Sj ^ lmaac _ £i maab + ^Ima bc - ^Rea ac - -f a d c ec) Rea ac - Im{P*(c}}, (33) 

d{a)/dt = -iCl/2(b) - jc/2(a) - iP, (34) 

d(b)/dt = i5 n {b) - iJ/2{&) - ifl/2(a) - j DX /2(b), (35) 

d(c)/dt = i(6 n - Sj)(c) - i J 12(b) - 1IX /2(c), (36) 



where ^ ab = 7c + Jdx , Ihc = Idx + lix and 7 ac = 7c + Jix arc reduced damping constants used in correlators. 

Additionally, we introduced terms corresponding to pure decoherence to the correlators, with rates being 7„j ec \ 
^(dec) ^(dec) ^ i-Qjjg a jj ows ^ acc0 unt temperature effects and other incoherent processes. 

The numerical calculation for dynamics of the modes ni X (t) and TiDx(t) calculated using system of Eqns. (3)-(5) 
written in the main text and advanced system (|28p ~ (|36[) presented above give the same result for the chosen pumping 
conditions, and thus ensures the validity of used mean field approximation. Therefore, while the model formulated 
in Section B is rather general, we do not use it for calculations and leave the problem of accounting of decoherent 
processes for future works. 



C: Simplified model 



Here we present model based on the simple three coupled oscillators Hamiltonian used in Ref. [8| . The Hamiltonian 
and corresponding wave function read 

/ ujix -J/2 \ /f/jf\ 
H = hi -J/2 Lj DX 0/2 , I*) = , (37) 

V 0/2 lo c J V *c / 

where the wave- function |\I>) of a dipolariton can be expanded in the basis of uncoupled states ^c, and &ix 

being field operators for cavity photon, direct exciton and indirect exciton, respectively. Here we essentially assume 
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FIG. 4: (Color online) (a), (b) and (c): Energy spectrum of dipolariton modes (solid lines) as a function of normalized electric 
field F/Fo plotted for Sn = —1.67. We assume equal IX-DX and C-DX coupling constants, measure energy in units of J and 
choose cavity mode as a reference point, Ec = 0. Pink dotted lines show energy position of the pumping of the cavity mode, (d) 
Frequency of long-term indirect exciton density oscillations v for different electric field, showing the possible tunable window of 
the signal frequency. Tunnelling frequency used for calculations is J/2n = 1.45 THz. (e) Amplitude of oscillations showing the 
correspondence between pumping conditions and mode detuning, (e) Quality factor of the oscillations defined as ratio between 
amplitude and decrement of oscillation. 



macroscopic occupation numbers for each mode, which allows us to use mean- field description and treat ^c, ^dx 
and ^ ix as classical fields. This enables describing the dynamics of the system by a closed set of time-dependent 
Schrodinger equations. To separate the rapidly oscillating part of dynamics we adopt the rotating frame using 
transformations: ^ IX {t) = ipix(t)e- iuJlxt , $ D x(t) = ipDx(t)e~ tUDxt and = ip c (t)e~ luJct Q. Finally, the 

system of equations describing dynamics of dipolariton system on a mean-field level reads 



dipix{t) J , ,., iS , t lix , u . 
— = t^4'Dx(t)e ldjt - — 1>i X (t), 



dipDx(t) 



- i^ix(t)e- i5jt - i^o(*)e- Wnt - i^^Dx(t), 



n 



Wait) 

at 



(38) 
(39) 
(40) 



where the photonic detuning 5q = loc ~ ^dx for the cavity photon coupled to the direct exciton, and energy distance 
between spatially indirect and direct excitons 5j = loj X — i^dx were introduced. The frequency of the resonant optical 
pumping Po(t) is detuned from the photon mode by A = up — ujq (second term in RHS of Eq. (|38p ). 

The presented description has shown to be well-suitable in numerous works devoted to physics of exciton-polaritons. 
However, it has serious drawback since does not account for decoherence in the real systems. The dissipation though 
can introduced in the system in phenomenological way, as we did adding last terms in Eqns. (|38j) . They are 
proportional to damping rates for each mode ji = 2ir/Ti, i = C, DX, IX. The lifetimes of the modes being typically 
Tc ~ 3 ps, t dx ~ 1 ns and tj X rj 100 ns. In the case where decoherence plays substantial role in the dipolariton 
system, this approach is not valid and one has to use the model from Section A based on the Master equation. 



10 




FIG. 5: (Color online) (a), (b) and (c): Energy spectrum of dipolariton modes (solid lines) as a function of normalized electric 
field F/Fo plotted for 8q = 0.5. Pink dotted lines show energy position of the pumping of the cavity mode, (d) Frequency 
of long-term indirect exciton density oscillations v for different electric field, showing the range where signal frequency can be 
tuned, (e) Amplitude of oscillations showing the correspondence between pumping conditions and mode detuning, (e) Quality 
factor of the oscillations defined as ratio between amplitude and decrement of oscillation. 



D: Dynamics analysis 

In the work we have studied the dynamics of indirect exciton density oscillations which induce emission of coherent 
THz radiation. These oscillations emerge due to the strong coupling between three modes, namely: indirect cxcitons 
(IX), direct excitons (DX) and cavity photons (C). In general, the dynamics of dipolaritons described by system of 
coupled equations for occupation numbers and correlators is complex due to the presence of mixing terms, detuning, 
decay and tunable pumping conditions. It can be described by characteristic energies HJ and Ml which denote IX-DX 
and C-DX coupling, respectively, detuning parameters Sj = u>ix — ^>dx and 6q = cue —^dx, and pumping frequency 
detuning from the cavity mode A = up — ujc- Consequently, in order to comprehend the dynamics of the system one 
needs to study the dependence of density oscillations on five independent parameters, which is a non-trivial task. 

The coupling constants of light-matter interaction Ml and indirect exciton to direct exciton coupling h J are intrinsic 
characteristics of the structure. The structure can be designed in such a way that these two characteristics coincide: 
hfl = HJ, in which case DX, IX and C modes arc efficiently intermixed. This regime has already been studied in the 
samples where dipolariton modes have been observed for the first time @, S] with Ml = HJ = 6 meV. Therefore, the 
behavior of dipolaritons can be controlled tuning the energy spacing between modes by altering the applied electric 
field (<5j) and tuning the incidence angle of the excitation light (Sq). 

We shall specifically consider two particular cases. In the first case the cavity mode is far detuned from the IX- 
DX resonance, which partially resembles the weak light- matter coupling (regime /). The second case refers to the 
resonance of all three modes achieved for a certain incidence angle where Sq is small (regime //) . For each regime we 
study the dependence of properties of indirect exciton density oscillations on the applied voltage. To keep the results 
general, the analysis is performed in dimensionless units F/Fq, where electric field is normalized to value where IX 
and DX resonances anti-cross, Fq. Moreover, we measure the frequency in units of tunnelling coupling J. Finally, the 
last parameter which governs the dynamics of the system is the pumping frequency ujp = A + wc . 

In the main text of this paper we present the plot of indirect indirect exciton density oscillations for two particular 
situations (Fig. 2 in the main text). In general case the dynamics exhibits turning-on transitional region lasting for 
several picoseconds after the arrival of the pulse followed by subsequent long-term oscillations of density decaying in 
time. Here we vary the pumping conditions and study the "quality" of oscillations. To do this we use the following 
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set of parameters. Complementary to the frequency of long-range oscillations, we find the averaged amplitude of the 
oscillations in transient regime AN to define the efficiency of pumping conditions. This quantity strongly affects the 
power output of generated initial THz signal and it is important for achievement of powerful pulsed THz emission. 
Additionally, we define dimensionless oscillation quality factor £ defined as a ratio of long-term oscillation amplitude 
AN to its decrement in time 6N. Therefore, high £ stands for slowly decaying oscillations with comparably large 
amplitude. 

We consider the first regime when cavity mode is far-detuned from the IX-DX anti-crossing region, Sq = —1.67. In 
energy diagrams in Fig. EJa,b,c) we chose three pumping energies corresponding to the resonant pumping (A = 1.92), 
blue detuned pumping (A = 2.5) and red-detuned pumping (A = —0.67). Solid lines there denote upper, middle 
and lower dipolariton branches (red, blue and green lines, respectively). The energy is counted from the bottom 
of the cavity mode (Ec = 0). Fig. [4(d) shows dependence of the frequency of mx oscillations as function of the 
electric field for different pumping energy shown in energy diagrams in Fig. Eta,b,c) (pink line). One can see that 
it has the minimum value for almost resonant electric field F '/ 'Fq ~ 0.9 and approximately corresponds to tunnelling 
coupling v sa J/2tt = 1.45 THz. Also, it increases with the increase of detuning between IX and DX modes, Sj. 
Moreover, we note that the frequency of stabilized oscillations does not depend on the pumping conditions. In Fig. 
HJc) we show amplitude of oscillations. It differs drastically for three pumping frequencies. One can see that AiV" 
increases if pump coincides with energy bare transition or dipolariton mode, and this correlation can be tracked in 
Fig. 0(c) for all detunings. However, calculations show that oscillations with the largest amplitude correspond to 
the strongest decrement of oscillations, which result into minimal £ in the resonant region (Fig. 2(f)). Out of the 
resonance both amplitude oscillation and decay are small, and the realistic region should be chosen between electric 
field corresponding to maxima of the quality factor £. 

Next, we study the second regime where cavity mode is detuned from direct exciton by 8q = 0.5. In the IX-DX 
resonance point dipolariton represents strong mixture of bare modes with two anti-crossings shown in Fig. |4(a-c). The 
dynamics of the system was calculated for three pumping energies which change from resonance (Fig. [1(a), A = 0) 
to out-of-resonance conditions (Fig. 0(c), A = —0.87). Fig. [5(d) shows no influence of pumping conditions on the 
oscillation frequency, while the position of the minimum shifts to F/Fq w 1.1 point. An amplitude of oscillations 
in Fig. [5(e) shows the same correlation with respect to pumping conditions observed in the first regime, while its 
absolute value is larger. However, the quality factor £ of oscillation in second regime is several times lower due to 
faster decay of oscillations (Fig. [3(f)). 

Thus, we can conclude that the first regime is appropriate for the cases where stable long-term emission of THz 
radiation is required, while the second regime is needed for the pulsed generation of powerful signals. 

E: Sketch of system with supplemental THz cavity 




MM (THz) MM (THz) 
DBR (o) 



FIG. 6: (Color online) Sketch of the geometry where power of THz emission is increased due to Purcell effect (top view). Here 
DBR (o) represent optical cavity mirrors (green), while MM (THz) denote metallic mirror cavity tuned to reflect signal of THz 
frequency (red). The red lines show peculiar polar pattern of dipolariton THz emitter. 
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